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ABSTRACT 



Identifying terrestrial planets in the habitable zones (HZs) of other stars is 
one of the primary goals of ongoing radial velocity and transit exoplanet surveys 
and proposed future space missions. Most current estimates of the boundaries 
of the HZ are based on 1-D, cloud-free, climate model calculations by Kasting 
et al. (1993). However, this model used band models which were based on older 
HITRAN and HITEMP line-by-line databases. The inner edge of the HZ in 
Kasting et al. (1993) model was determined by loss of water, and the outer edge 
was determined by the maximum greenhouse provided by a CO2 atmosphere. 
A conservative estimate for the width of the HZ from this model in our Solar 
system is 0.95-1.67 AU. 

Here, an updated 1-D radiative-convective, cloud-free climate model is used 
to obtain new estimates for HZ widths around F, G, K and M stars. New H2O 
and CO2 absorption coefficients, derived from the HITRAN 2008 and HITEMP 
2010 line-by-line databases, are important improvements to the climate model. 
According to the new model, the water loss (inner HZ) and maximum greenhouse 
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(outer HZ) limits for our Solar System are at 0.99 AU and 1.70 AU, respectively, 
suggesting that the present Earth lies near the inner edge. Additional calculations 
are performed for stars with effective temperatures between 2600 K and 7200 K, 
and the results are presented in parametric form, making them easy to apply 
to actual stars. The new model indicates that, near the inner edge of the HZ, 
there is no clear distinction between runaway greenhouse and water loss limits 
for stars with Tg// < 5000 K which has implications for ongoing planet searches 
around K and M stars. To assess the potential habitability of extrasolar terrestrial 
planets, we propose using stellar flux incident on a planet rather than equilibrium 
temperature. This removes the dependence on planetary (Bond) albedo, which 
varies depending upon the host star's spectral type. We suggest that conservative 
estimates of the HZ (water loss and maximum greenhouse limits) should be used 
for current RV surveys and Kepler mission to obtain a lower limit on t]^, so that 
future flagship missions like TPF-C and Darwin are not undersized. Our model 
does not include the radiative effects of clouds; thus, the actual HZ boundaries 
may extend further in both directions than the estimates just given. 

Subject headings: stars: planetary systems 

1. Introduction 

As of November, 2012, more than 800 extra-solar planetary systems have been detected^, 
and > 2000 additional candidate systems from the i^e|?/er mission are waiting to be confirmed 
(Batalha et al. 2012). One of the primary goals of the ongoing radial velocity (RV) and transit 
surveys is to identify a terrestrial mass planet (0.3 — lOM®) in the so-called Habitable Zone 
(HZ), which is traditionally defined as the circumstellar region in which a terrestrial- mass 
planet with a CO2-H2O-N2 atmosphere can sustain liquid water on its surface^ (Huang 1959; 
Hart 1978; Kasting et al. 1993; Underwood et al. 2003; Selsis et al. 2007b; Kaltenegger et al. 
2011b; Pierrehumbert & Gaidos 2011). Several potential HZ planet candidates have already 
been detected, (Udry et al. 2007; Pepe et al. 2011a; Borucki et al. 2011; Bonfils et al. 2011; 
Borucki et al. 2012; Vogt et al. 2012; Tuomi et al. 2012a,b) and it is expected that this 
number will greatly increase as time passes (Batalha et al. 2012). In the near future we may 
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^Abe et al. (2011) studied habitability of water-limited 'land' planets and found that they could remain 
habitable much closer to their host stars. However, Abbot et al. (2012) found that a waterworld would have 
a narrower HZ owing to lack of weathering-climate feedback. 
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be able to study habitable planets orbiting nearby M stars. These planets are relatively close 
to their parent stars, leading to shorter orbital periods and an increase in the probability 
of a transit. NASA's James Webb Space Telescope (JWST), scheduled to launch in 2018, is 
considered to be marginally capable of obtaining a transit spectrum of an Earth-like planet 
orbiting a late M dwarf (Clampin et al. 2007; Kaltengger & Traub 2009; Deming et al. 2009). 
Several other surveys are either underway (Nutzman & Charbonneau 2008, MEARTH) or 
getting ready to be commissioned (Mahadevan et al. 2012, HPF) in an attempt to discover 
rocky planets in the HZs of low mass stars. 

The HZ limits that were cited in many recent discoveries were obtained from 1-D 
radiative-convective, cloud-free climate model calculations by Kasting et al. (1993). For 
our Sun, these authors estimated the boundaries of the HZ to be 0.95 AU for the inner edge 
and 1.67 AU for the outer edge. These values represent the "water loss" and "maximum 
greenhouse" limits, respectively. Other, less conservative limits for the inner edge are the 
"runaway greenhouse" and "recent Venus" limits. The latter estimate is empirical, based on 
the inference that Venus has not had liquid water on its surface for at least the last 1 billion 
years (Solomon & Head 1991). For the outer edge, there is a corresponding "early Mars" 
empirical estimate, based on the inference that Mars did have liquid water on its surface 3.8 
billion years ago. (The "1st CO2 condensation" limit of Kasting et al. (1993), should now be 
disregarded, as it has been shown that CO2 clouds generally warm a planet's climate (Forget 
& Pierrehumbert 1997)). Some studies have investigated the effects of clouds on planetary 
emission spectra of Earth-like planets in a ID model (Kitzmann et al. 2011a,b), while others 
studied the habitability of specific systems, particularly Gl 581, in ID (Wordsworth et al. 
2010; Von Paris et al. 2011a; Kaltenegger et al. 2011a) and 3D (Wordsworth et al. 2011; 
Pierrehumbert 2011). Several other studies (Underwood et al. 2003; Selsis et al. 2007b) 
parameterized these results to estimate relationships between HZ boundaries and stellar 
parameters for stars of different spectral types. 

Although these studies provided useful estimates of the HZ width, the Kasting et al. 
(1993) model has become outdated, for several reasons: 

1. Kasting et al. (1993) used 'band models'^ for H2O and CO2 absorption in the thermal- 
infrared. These coefficients were considered valid up to ~ 700 K. These coefficients 
were later replaced (Mischna et al. 2000) by coefficients generated using the correlated-/;; 
technique (Mlawer et al. 1997; Kato et al. 1999). A line-by-line (LBL) radiative transfer 
model, in this case LBLRTM (Clough & locono 1995), was used to generate detailed 
spectra for H2O and CO2 at a variety of different temperatures and pressures. Once 



■^See Appendix B of Kasting (1988) for a detailed description of the band model. 



-4- 



the detailed spectra were calculated, separate broad-band k-coefficients for both H2O 
and CO2 were generated by R. Freedman using standard procedures. But these coeffi- 
cients were only derived for temperatures < 350 K and should therefore underestimate 
thermal-IR absorption in warm, moist greenhouse atmospheres. (This prediction was 
verified by direct experimentation with that model.) Furthermore, the coefficients 
adopted by Mischna et al. (2000) and used in subsequent climate modehng studies by 
the Kasting research group were obtained using HITRAN 1996 database and had not 
been updated since then. 

2. Recent studies (Halevy et al. 2009; Wordsworth et al. 2010) have pointed that the 
Kasting et al. (1993) model may have significantly overestimated absorption of thermal- 
IR radiation by collision-induced absorption (CIA) bands of CO2, which may affect the 
outer edge of the HZ. 

3. The Kasting et al. (1993) calculations spanned stellar effective temperatures from 7200 
K to 3700 K, corresponding approximately to stellar classes FO to MO. Stellar effective 
temperature affects the HZ boundaries because the radiation from F stars is bluer rela- 
tive to that from the Sun, whereas the radiation from K and M stars is redder, and this 
affects calculated planetary albedos. The HZ limits from Kasting et al. (1993) model 
do not include M stars with effective temperatures lower than 3700 K. As pointed out 
above, such stars are promising candidates for current observational surveys because 
their HZs are closer to the star. Therefore, potential rocky planets in the HZs will have 
shorter orbital periods and higher probability of transit. 

In this paper we address all the above major issues with the goal of deriving new, 
improved estimates for the boundaries of the HZ. The outline of the paper is as follows: In 
§2 we describe our 1-D cloud-free climate model, corresponding model updates and model 
validation with other studies. In §3 we present results from our climate model and discuss 
various HZ limits for our Earth. §4 presents HZ boundaries around F, G, K and M spectral 
stellar spectral types, then provides a generalized expression to calculate HZ boundaries and 
compares these boundaries with previous studies. We discuss the implications of these new 
results for currently known exoplanet planetary systems in §5 and present our conclusions 
in §6. 

2. Model description 

We used a one-dimensional, radiative-convective, cloud-free climate model based on 
Kasting (1988) for the inner edge of the HZ (IHZ) and Kasting (1991) for the outer edge 



- 5 - 



of the HZ (OHZ) calculations. Following Kasting et al. (1993), we assumed an Earth- 
mass planet with an H2O (IHZ) or CO2 (OHZ) dominated atmosphere for our base model. 
Sensitivity studies for different planetary masses are described in the following section. Both 
the inner and outer edge calculations relied on so-called "inverse climate modeling" , in which 
the surface temperature is specified, and the model is used to calculate the corresponding 
solar flux needed to sustain it. To do this, the atmosphere was divided into 101 layers, and 
a specific pressure-temperature profile was assumed. For the inner edge, this consisted of a 
moist pseudoadiabat extending from the surface up to an isothermal (200 K) stratosphere. 
Methodology for calculating the pseudoadiabat was taken from Appendix A of Kasting 
(1988). The surface temperature was varied from 200-2200 K during the course of the 
calculations. For the outer edge, the surface temperature was fixed at 273 K, and the CO2 
partial pressure was varied from 1 to 37.8 bar (the saturation CO2 partial pressure at that 
temperature). A moist H2O adiabat was assumed in the lower troposphere, and a moist CO2 
adiabat was used in the upper troposphere when condensation was encountered, following 
the methodology in Appendix B of Kasting (1991). 

H2O and CO2 clouds are neglected in the model, but the effect of the former is ac- 
counted for by increasing the surface albedo, as done in previous climate simulations by the 
Kasting research group (Kasting 1991; Haqq-Misra et al. 2008). It has been argued that this 
methodology tends to overestimate the greenhouse effect of dense CO2 atmospheres (Gold- 
blatt & Zahnle 2011). By contrast, our neglect of CO2 clouds may cause us to underestimate 
the greenhouse effect of such atmospheres (Forget & Pierrehumbert 1997). Realistically de- 
termining the effects of clouds would require a 3-D climate model, as most clouds form in 
updrafts, which are absent in 1-D models. Some 1-D climate modeling studies include par- 
tial cloud coverage (Selsis et al. 2007b) and/or parameterized microphysical cloud model 
(Colaprete & Toon 2003; Zsom et al. 2012), but we do not consider them here because we 
can not model them self-consistently in our model. The effects of clouds on the inner and 
outer edge boundaries are qualitatively understood, as discussed later in the paper. Testing 
these predictions quantitatively using 3-D climate models should be a fruitful topic for future 
research. 

Radiative transfer was handled by methods used in recent versions of the Kasting group 
climate model but with updated absorption coefficients (see next section). A 6 two-stream 
approximation (Toon et al. 1989) was used to calculate the net absorbed solar radiation 
for each of the 101 layers, using separate eight-term, correlated-A; coefficients for both CO2 
and H2O to parameterize absorption in each of the 38 solar spectral intervals ranging from 
0.2 — 4.5/im. These terms are convolved with each other in each spectral interval, resulting 
in 64 separate radiative transfer calculations per interval. The solar flux was averaged over 
six zenith angles (11.0°, 25.3°, 39.6°, 54°, 68.4°, 82.8°) using Gaussian quadrature. The net 
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outgoing infrared radiation per layer was calculated using separate eight-term correlated- /c 
coefficients for H2O and CO2 in 55 spectral intervals extending from — 15,000 cm~^. We 
used double gauss quadrature in place of a standard gaussian scheme (Sykes 1952; Thomas 
& Stamnes 2002) using a code written by Ramirez. Half of the k-coefficients are chosen 
within the g-space interval 0.95-1.00 for improved resolution of the steeply rising portion of 
the cumulative distribution function, yielding smoother stratospheric temperature behavior. 

These coefficients also needed to be convolved with each other, as in the solar calculation. 
This produces 8x 8 x 55 = 3520 separate thermal-IR radiative transfer calculations at each 
time step in the climate model. This number is multiplied by a factor of 6 when we include 
CH4 in the model, using 6-term sums, and by another factor of 6 when we include C2H6. 
Thus, from a practical standpoint, the utility of this approach diminishes as the number of 
included greenhouse gases increases. 

2.1. Model Updates 

The following are the most significant updates to the climate model: 

1. We have derived new k- coefficients using a tool called KSPECTRUM. It is a program 
to produce high-resolution spectrum of any gas mixture, in any thermodynamical con- 
ditions, from line-by-line (LBL) databases such as HITRAN 2008 (Rothman et al. 
2009) and HITEMP 2010^ (Rothman et al. 2010). It is intended to produce rehable 
spectra, which can then be used to compute k-distribution data sets that may be used 
for subsequent radiative transfer analysis. The source code and a detailed description 
of the program is available at http://code.google.eom/p/kspectruin/. 

We have produced two sets of coefficients, one using HITRAN 2008 and another us- 
ing the HITEMP 2010 database. For the HITRAN database we generated a matrix of 
8-term absorption coefficients for both H2O and CO2, using KSPECTRUM, for the fol- 
lowing range of pressures and temperatures: p{bar) = [10~^, 10~^, 10~^, 10~^, 10~^, 1, 10, 10^] 
and T{K) = [100,150,200,250,300,350,400,600]. In the case of HITEMP, 8-term 
absorption coefficients were derived only for H2O, as our IHZ is H20-dominated at 
high temperatures (> 300 K) with only trace amounts of CO2 (330 parts per mil- 
lion). The following grid was used to derive the H2O HITEMP coefficients: p{bar) = 
[10-\1,10,102] and T{K) = [350,400,600]. The grid is condensed because of the 
high number of line transitions in the HITEMP database compared to HITRAN. The 
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computational resources needed to derive absorption coefficients for the entire range of 
pressure & temperatures would be prohibitively large. Moreover, as discussed further 
below, we justify the selection of this condensed grid by showing that the differences 
in coefficients generated from HITRAN and HITEMP become negligible below 350 K. 

In generating the /c— coefficients, we have used different methodologies for CO2 and 
H2O. For C02, we truncated the spectral lines at 500 cm~^ from the line center. Ex- 
perimental evidence indicates that the absorption by CO2 is overestimated if Lorentzian 
line shapes are used (Burch et al. 1969; Fukabori et al. 1986; Bezard et al. 1990; Halevy 
et al. 2009). Therefore, we used the prescription of Perrin & Hartmann (1989) for 'sub- 
Lorentzian' absorption in the far wings of the lines when running KSPECTRUM. For 
H2O, we truncated the spectral lines 25 cm~^ and overlaid a semi-empirical "contin- 
uum absorption". The Lorentz line shape is known to underestimate absorption for 
H2O in the far wings (Halevy et al. 2009), possibly because of the tendency of H2O to 
form dimers. The corresponding continuum absorption is therefore 'super-Lorenztian' 
for H2O, and we have used the 'BPS' formalism of Paynter & Ramaswamy (2011) to 
parameterize this absorption. 

2. We have included Rayleigh scattering by water vapor, as it can become important for 
wavelengths up to l/im (which is where the Wien peak occurs for low mass stars). 
Rayleigh scattering by water was also considered by Kasting (1988); Kasting et al. 
(1993), but these authors used the scattering coefficient for air because the coefficient 
for H2O was not available, or at least not known to them. The following expression 
for the scattering crosssection was adopted (Allen 1976; Vardavas & Carver 1984; Von 
Paris et al. 2010): 

aK,H.o(A) = 4.577 X 10-- (^±1^) (J) cm^ (1) 

Here, D is the depolarization ratio (0.17 for H2O, Marshall & Smith (1990)), r is 
the wavelength (A)-dependent refractivity which is calculated as r = 0.85rdryair (Edlen 
1996), Tdryair IS obtained from Eq.(4) of Bucholtz (1995), and A is in microns. By 
comparison, Selsis et al. (2007a) used a H2O Rayleigh scattering cross-section of 2.32 x 
10"^'' cm^ at 0.6 /im. Evaluating Eq.(l) at 0.6 yum gives a value of 2.6 x 10"^^ cm^, 
which is similar to the Selsis et al. (2007a) value. 

3. Previous climate model calculations by our group and others (Kasting et al. 1984; 
Pollack et al. 1987; Kasting 1991; Forget & Pierrehumbert 1997; Mischna et al. 2000) 
parametrized collision-induced absorption (CIA) by CO2 near 7 /im and beyond 20 /im 
by the formulation given in the Appendix of Kasting et al. (1984). This process is an 
important source of thermal-IR opacity in the types of dense, C02-rich atmospheres 
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predicted to be found near the outer edge of the habitable zone. In our model we have 
updated CO2-CIA using the parametrization described in Gruszka & Borysow (1997); 
Baranov et al. (2004); Halevy et al. (2009). 

4. The Shomate Equation^ was used to calculate new heat capacity (cp) relationships for 
CO2 and H2O. Notably, at low temperatures, the heat capacity for CO2 decreased by 
~ 30% relative to values in our previous model. This increased the dry adiabatic lapse 
rate, g/cp where g is gravity, by an equivalent amount but had surprisingly little effect 
on computed surface temperatures, apparently because the steeper lapse rate in the 
upper troposphere was largely compensated by a decrease in tropopause height. See 
Ramirez et al. (2012a) for further details. 

2.2. Model Validation 

We have checked the accuracy of our climate model by comparing the output both with 
published results and with the 1-D line- by-line radiative transfer model SMART {Spectral 
Mapping Atmospheric Radiative Transfer) developed by D. Crisp (Meadows & Crisp 1996; 
Crisp 1997). SMART is a well-tested model (Robinson et al. 2011) which accesses some of 
the same databases as does KSPECTRUM; however, its development and implementation 
are entirely independent. By comparing specific cases of interest with SMART, we can 
gain confidence that our calculated fluxes are correct, or at least that they are consistent 
with our assumptions about CO2 and H2O line shapes. For all our climate models that 
are compared with SMART, we used 70 atmospheric layers (we use 101 layers for all our 
HZ calculations). We could not use 101 layers in our flux comparisons due to numerical 
accuracy issues with SMART at high enough vertical resolution, although 70 layers produced 
a sufficiently accurate result with SMART. 

2.2.1. Dense CO2 atmosphere 

Dense C02-rich atmospheres have been suggested as warming agents for early Mars 
(Pollack et al. 1987; Kasting 1991; Forget & Pierrehumbert 1997; Tian et al. 2010). Planets 
close to the outer edge of the HZ may develop dense, C02-rich atmospheres as a consequence 
of outgassing from volcanism, which can only be balanced by surface weathering if the planets 
surface temperature remains above freezing. The CO2 feedback effect fails at some distance 



■■'http : //webbook . nist . gov/cgi/cbook . cgi?ID=C124389&Unit s=SI&Mask=l#Thermo-Gas 



because CO2 begins to condense out of the atmosphere, lowering the tropospheric lapse 
rate and reducing the greenhouse effect. CO2 is also an effective Rayleigh scatterer (2.5 
times better than air), and so a dense CO2 atmosphere is predicted to have a high albedo, 
which offsets its greenhouse effect (Kasting 1991). The OHZ boundary can then be taken as 
this "maximum greenhouse limit" where Rayleigh scattering by CO2 begins to outweigh the 
greenhouse effect. 

Fig. 1 shows net outgoing long-wave radiation (OLR) versus wavenumber in the range 
0—2000 cm~^ for a Mars- mass planet with a 2-bar CO2 atmosphere and a surface temperature 
of 250 K. The solar constant is assumed to be 75% of its present value (1360 Wm~^), 
matching the solar flux incident on early Mars (3.8 Gyr). The integrated flux over all bands 
at the top of the atmosphere from our model (86 Wm~^, blue solid curve) matches well 
with SMART (88.4 Wm~^, dashed red curve). Our model has a coarser spectral resolution 
than does SMART, and it appears that between 800 and 1200 cm~^ the differences could 
become important. But this is compensated by the fact that our OLR in these intervals 
can be considered as a running average of the OLR from SMART. Nevertheless, most of 
the difference in the OLR arises from significant absorption in the 667 cm~^ (15 micron) 
vibrational band of CO2 which is closer to the peak of the black-body curve. 

A similar study for early Mars conditions with a 2-bar CO2 atmosphere was considered 
by Wordsworth et al. (2010). They have also used KSPECTRUM to derive their absorption 
coefficients and truncated the spectral lines at 500 cm~^ from the line center for CO2, as 
done here. As our surface albedo for this calculation (0.2) is also the same, we can directly 
compare the results from both studies. Fig. 2c from Wordsworth et al. (2010) shows that 
the net OLR from their model is 88.17 Wm~^ compared to our 86 Wm~^. The differences 
are due to the different number of atmospheric layers used in these models. Wordsworth et 
al. (2010) used 22 layers in their model, compared to 70 layers in our SMART comparison 
climate models. The number of vertical atmospheric layers used in the model affects the 
OLR because the Toon et al. (1989) algorithm, used in both the models, assumes that each 
layer is isothermal. With few isothermal layers, more IR radiation is emitted from the upper 
part of each layer, which is a little hotter than it should be and which has the smallest optical 
depth, as measured from the top of the atmosphere^. 



'We ran our climate model with 22 layers and found that our OLR increased to 89.1 Wm 
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Fig. 1. — Plot of outgoing long- wave radiation vs. wavenumber for the — 2000cm~^ region 
comparing our OLR (blue solid curve) to that from SMART (red dashed curve). This 
calculation is for early Mars conditions, 2 bar CO2, and constant stratospheric and surface 
temperatures of 167 and 250 K, respectively. The corresponding 250 K black body curve is 
shown in black. The integrated flux over all bands at the top of the atmosphere is 86 Wm~^ 
for our model and 88.4 Wm-^ for SMART. 



2.2.2. Dense H2O atmosphere 



The inner edge of the HZ in our model is determined by the so-called "moist green- 
house effect", in which the stratosphere becomes water-dominated, leading to rapid escape 
of hydrogen to space. Fig. 2 shows the net outgoing IR as a function of wavenumber for a 
dense H2O atmosphere. Here, we assumed an Earth-mass planet with a surface temperature 
of 400 K and a surface albedo of 0.3. The stratospheric temperature is assumed to be con- 
stant at 200 K. The stratosphere becomes tenuous at these high surface temperatures and 
has little effect on the outgoing IR flux. The background gas is 4 bar of N2 and the total 
surface pressure is 6.5 bar (These conditions were assumed for specific intercomparison with 
SMART for this test case). The flux incident at the top of the atmosphere is assumed to be 
the current solar flux at Earth's distance from the Sun. 

As with the dense CO2 case, in Fig. 2 we compare our model (solid blue curve) with 
SMART (dashed red curve) for the dense H2O atmosphere. Although both model spectra ap- 
pear to be in good agreement, the integrated flux over all bands at the top of the atmosphere 
from our model is 285 Wm~^ compared to 297 Wm~^ from SMART. The differences arise in 
the window region of the water vapor (800 — 1200 cm~^) and also in between 300 — 600 cm~^, 
where our model absorbs more than SMART. A possible reason is that we are using the BPS 
continuum, as opposed to the 'CKD' continuum (Clough et al. 1989) used by SMART. The 
BPS formalism is based on empirical measurements which take into account the contribution 
of dimers, resulting in more absorption of outgoing IR radiation (see Paynter & Ramaswamy 
(2011) Table 3; Shine et al. (2012)). 

3. Results 

In the subsections that follow, we estimate HZ boundaries around a star similar to our 
Sun. We first compare results from our model using HITRAN and HITEMP databases, 
estimate HZ limits for non-Earth-like planets and discuss the effect of clouds on the HZ 
boundaries. 



3.1. Inner Edge of the HZ (IHZ) 

The inner edge of the HZ is calculated by increasing the surface temperature of a fully 
saturated "Earth" model from 220 K up to 2200 K. The effective solar flux S^ff, which is the 
value of solar constant required to maintain a given surface temperature, is calculated from 
the ratio between the net outgoing IR flux F/r and the net incident solar flux Fsol, both 
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Fig. 2. — Plot of outgoing long-wave radiation vs. wavenumber for the — 2000 cm~^ 
region comparing our OLR (blue solid curve) to that from SMART (red dashed curve). This 
calculation is for a dense H2O atmosphere with stratospheric and surface temperatures of 
200 K and 400 K, respectively. The corresponding 400 K blackbody curve is shown in black. 
Both the models appear to be in good agreement. The integrated flux over all bands at 
the top of the atmosphere is 285 Wm~^ for our model and 297 Wm~^ for SMART, with 
differences possibly arising from different formalisms of continuum absorption (BPS versus 
CKD. See text for details). 
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evaluated at the top of the atmosphere. The total flux incident at the top of the atmosphere 
is taken to be the present solar constant at Earth's orbit 1360 Wm~^. The planetary albedo 
is calculated as the ratio between the upward and downward solar fluxes. 

The calculated radiative fluxes, planetary albedo and water vapor proflle for various 
surface temperatures are shown in Fig. 3. Absorption coefficients derived from the HITEMP 
2010 database, overlaid by BPS formalism (Paynter & Ramaswamy 2011), were used in 
generating these results. Fig. 3(a) shows that -F/^j increases with surface temperature and 
then levels out at 291 Wm~^, as the atmosphere becomes opaque to infrared radiation at all 
wavelengths^. Beyond 2000 K, Fm increases again as the lower atmosphere and surface begin 
to radiate in the visible and near-IR, where the water vapor opacity is low. Fsol initially 
increases as a consequence of absorption of near-IR solar radiation by H2O. It then decreases 
to a constant value (264 Wm~^) at higher temperatures as Rayleigh scattering becomes 
important. Planetary albedo (Fig. 3(b)) provides an alternative way of understanding this 
behavior. It goes through a minimum at a surface temperature of 400 K, corresponding to 
the maximum in Fsol-, and then flattens out at a value of 0.193. 

The inner edge of the HZ for our Sun can be calculated from Fig. 3(c). The behavior 
of Fir and Fsol causes Sejf to increase initially and then remain constant at higher tem- 
peratures. Two limits for the IHZ boundary can be calculated. The flrst one is the "moist 
greenhouse" (or water-loss) limit which is encountered at a surface temperature of 340 K 
when Seff = 1.015. At this limit, the water vapor content in the stratosphere increases 
dramatically, by more than an order of magnitude, as shown in Fig. 3(d). This is the rele- 
vant IHZ boundary for habitability considerations, although it should be remembered that 
the actual inner edge may be closer to the Sun if cloud feedback tends to cool the planet's 
surface, as expected.^ The orbital distance corresponding to the cloud-free water loss limit 
is d = 1/ S^fj = 0.99 AU for an Earth-like planet orbiting the Sun. 

The second IHZ limit is the runaway greenhouse at which the oceans evaporate entirely. 
The limiting Seff from Fig. 3(c) is 1.06 which corresponds to a distance of 0.97 AU. Both 
calculated IHZ limits are signiflcantly farther from the Sun than the values found by Kasting 
et al. (1993) (0.95 AU for the water-loss limit and 0.84 AU for the runaway greenhouse). 



^This value of 291 Wni ^ closely matches with the value from Fig. 4.37 of Pierrehumbert (2010) for a 
planet saturated with pure water vapor atmosphere and with a surface gravity of lOms"^. 

^The total H2O inventory assumed here is equal to the amount of water in Earth's oceans - 1.4 x 10^'* 
grams. This amounts to 2 x 10^^ atoms per cm~^. Once the stratosphere becomes wet, water vapor photolysis 
releases hydrogen which can escape to space by diffusion limited escape rate. The time scale for water loss 
approaches the age of the Earth when the mixing ratio of water is ^ 3 x 10"'^, which happens at a surface 
temperature of 340 K. 



The difference is caused by increased atmospheric absorption of incoming solar radiation by 
H2O in the new model. As pointed out by Kasting et al. (1993), a third estimate for the 
IHZ boundary can be obtained from radar observations of Venus by Magellan spacecraft, 
which suggest that liquid water has been absent from the surface of Venus for at least 
1 Gyr (Solomon & Head 1991). The Sun at that time was ~ 92% of the present day 
luminosity, according to standard stellar evolutionary models (Baraffe et al. 1998; Bahcall 
et al. 2001, See Table 2). The current solar flux at Venus distance is 1.92 times that of 
Earth. Therefore, the solar flux received by Venus at that time was 0.92 x 1.92 = 1.76 times 
that of Earth. This empirical estimate of the IHZ edge corresponds to an orbital distance 
of d = (1/1.76)°'^ = 0.75 AU for the present day. Note that this distance is greater than 
Venus' orbital distance of 0.72 AU because the constraint of surface water was imposed at 
an earlier time in the planet's history. 

3.2. Comparison of inner edge results using the HITEMP and HITRAN 

databases 

In Fig. 4 we show Fm as a function of surface temperature (similar to Fig. 3(a)). We 
wish to compare the outgoing IR calculated from HITRAN & HITEMP databases with and 
without overlaying the continuum absorption. Fig. 4 shows two significant differences: 

1. The limiting value of Fj^ which leads to a runaway greenhouse happens at a much 
higher value (440 Wm~^, black & green curves) when the BPS H2O continuum formal- 
ism is not implemented, and at a lower Fjr (291 Wm~^, red & blue curves) when the 
BPS continuum is included in our model. The continuum is based on measurements 
of absorption in the water vapor window regions (i.e 800 — 1200 cm~^ and 2000 — 3000 
cm~^). At high temperatures, the contribution of the continuum absorption in these 
window regions becomes significant, and this, in turn, decreases the outgoing IR flux. 

2. The moist-greenhouse (water loss) limit moves much closer to the Sun (to 0.87 AU) 
when continuum absorption is not included, as compared to 0.99 AU when it is included 
in our model. This is a direct consequence of the differences in Fm described above: 
When Fir increases with the continuum turned off, S^ff (ratio of Fjr to Fsol) increases 
and the IHZ distance d = l/S^f^ decreases. The result can be understood physically 
by noting that in the model where the continuum is absent, the planet needs more 
effective solar flux to maintain a given surface temperature because more thermal-IR 
radiation leaks away into space; hence, the IHZ boundary must move inward. 

A similar change can be seen in the runaway greenhouse limit: the 'No BPS' model 
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Fig. 3. — Inner edge of the habitable zone calculations from our updated climate model. 
Various parameters are shown as a function of surface temperature: (a) Net outgoing IR flux 
and net incident solar flux (b) planetary albedo (c) Effective solar flux, Sejf = Fj^/Fsol, 
and (d) water vapor profile. These calculations were performed with the HITEMP 2010 
database. The water-loss (moist greenhouse) limit, which is most relevant to habitability, 
is at 0.99 AU and runaway greenhouse is at 0.97 AU. The corresponding estimates from 
Kasting et al. (1993) climate model are 0.95 AU for the moist greenhouse and 0.84 AU for 
the runaway greenhouse. 
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transitions to runaway at a higher Sejf than does the 'with BPS' modeL The corre- 
sponding runaway greenhouse hmit changes from 0.97 AU (with continuum absorption) 
to 0.76 AU (without continuum absorption). Fig. 4 also shows that the upturn in Fjn 
beyond 800 K happens at lower surface temperatures when continuum absorption is 
not included. It should be remembered that this upturn in Fir happens because, as 
the surface warms, the region in the troposphere over which the temperature profile 
follows a dry adiabat expands upward, while the moist convective layer in the upper 
troposphere becomes thinner. Eventually, when the moist convective region (the cloud 
layer) begins thin enough, radiation emitted from the dry adiabatic portion of the at- 
mosphere begins to escape to space. The dry adiabatic lapse rate is steeper than the 
moist adiabatic lapse rate by about a factor of 9 (~ 10 K/km vs. 1.1 K/km); hence, 
the emitted radiation flux is much higher. This can be understood from the integrated 
form of Schwarzchild's equation, which shows that the emitted flux is proportional to 
the temperature gradient. (See, e.g., eq. A4 in Kasting et al. (1984)) Unlike Kasting 
(1988), we find that the emitted flux increases at all thermal-IR wavelengths shorter 
than 4 fim. The amount of visible radiation emitted remains negligible for surface tem- 
peratures of 2200 K or below. Without the continuum there are fewer lines to cause 
absorption in these thermal-IR bands and hence a lower temperature would suffice to 
cause the upturn. Fig. 4 also shows that the model that includes both HITEMP & 
continuum (red curve) is the one that absorbs the most outgoing IR radiation (which 
is the one that was used to derive inner HZ limits in §3.1). 



3.3. Outer Edge of the HZ (OHZ) 

In determining the outer edge of the HZ, the surface temperature of an Earth- like 
planet with 1-bar N2 atmosphere was fixed at 273 K and the atmospheric CO2 partial 
pressure, PCO2, was varied from 1 to 35 bars (the saturation vapor pressure for CO2 at that 
temperature). The stratospheric temperature was chosen as follows: The model atmosphere 
(Mars-like planet) in which the onset of CO2 condensation occurs has a cold-trap temperature 
of 154 K at an altitude where the ratio of the saturation vapor pressure to the ambient 
pressure is unity. We replace the temperature profile above this altitude with a constant 
temperature of 154 K. This allows us to calculate the solar fiux {Se/f) required to maintain 
a global mean surface temperature of 273 K as explained in §3.1. Our working hypothesis 
is that atmospheric CO2 would accumulate as these planets cooled because of the negative 
feedback provided by the carbonate-silicate cycle. Results from our model calculations are 
shown in Fig. 5. 
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Fig. 4. — Comparison of outgoing IR radiation (-F/zj) from HITEMP and HITRAN databases, 
with (blue & red curves) and without (black & green curves) continuum absorption. F/^ 
is lower in the 'with BPS' case for both the databases because the continuum absorption 
becomes significant in the water vapor window regions. 



The incident solar (Fsol) and outgoing IR (-F/_r) fluxes are shown in Fig. 5(a). Fjn 
decreases initially as CO2 partial pressure is increased; this is an indication of greenhouse 
effect of CO2. At ~ 10 bars, Fm asymptotically approaches a constant value as the atmo- 
sphere becomes optically thick at all infrared wavelengths. Fsol decreases monotonically 
with increases in CO2 partial pressure as a result of increased Rayleigh scattering. Corre- 
spondingly, the planetary albedo increases to high values at large CO2 partial pressures, as 
shown in Fig. 5(b). The solar and IR fluxes, acting in opposite directions, create a mini- 
mum of Seff = 0.325 at a CO2 partial pressure of ~ 8 bar (Fig. 5(c)), corresponding to a 
distance d = 1.70 AU. This defines the maximum greenhouse limit on the outer edge of the 
HZ. By comparison, Kasting et al. (1993) model predicted d = 1.67 AU for the maximum 
greenhouse limit. As emphasized earlier, radiative warming by CO2 clouds is neglected here, 
even though they should be present in this calculation. Therefore, our OHZ limit should be 
considered as a conservative estimate, that is, the real outer edge is probably farther out. 

As with the inner edge model, a more optimistic empirical limit on the OHZ can be 
estimated based on the observation that early Mars was warm enough for liquid water to 
flow on its surface (Pollack et al. 1987; Bibring et al. 2006). Assuming the dried up riverbeds 
and valley networks on martian surface are 3.8 Gyr old, the solar luminosity at that time 
would have been ~ 75% of the present value (See Eq.(l) in Gough (1981) and Table 2 in 
Bahcall et al. (2001)). The present-day solar flux at Mars distance is 0.43 times that of 
Earth. Therefore, the solar flux received by Mars at 3.8 Gyr was 0.75 x 0.43 = 0.32 times 
that of Earth. The corresponding OHZ limit today, then, would he d = (1/0.32)°'^ ^ 1.77 
AU. 

Note that this distance exceeds the maximum greenhouse limit of 1.70 AU estimated 
above indicating that to keep early Mars wet, additional greenhouse gases other CO2 and H2O 
may be required. In fact, Ramirez et al. (2012a) show that a 3-bar atmosphere containing 
90 percent CO2 and 10 percent H2 could have raised the mean surface temperature of early 
Mars above the freezing point of water. The warming is caused by the collision-induced 
absorption due to foreign-broadening by molecular hydrogen. It should be acknowledged 
that some authors (e.g., Segura et al. (2002, 2008)) do not agree that early Mars must have 
been warm; however, in our view, these cold early Mars models do not produce enough 
rainfall to explain valley formation (Ramirez et al. 2012a). 

3.4. Effect of clouds on the HZ boundaries 

We summarize various cloud-free HZ boundary estimates for Earth in Table 1. Although 
we updated our radiative transfer model to incorporate new absorption coefficients, this by 
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Fig. 5. — Outer edge of the habitable zone calculations from our climate model, shown as 
a function of CO2 partial pressure PCO2: (a) Net outgoing IR flux and net incident solar 
flux (b) planetary albedo (c) Effective solar flux. The maximum greenhouse limit, where the 
atmosphere becomes opaque to outgoing IR radiation, is at 1.70 AU (5'e// = 0.343). The 
previous estimate from Kasting et al. (1993) climate model was 1.67 AU. 
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itself may not yield a significantly better estimate for the width of the habitable zone. The 
reason is that it is widely acknowledged that the HZ boundaries will be strongly influenced 
by the presence of clouds. H2O clouds should move the inner edge inwards (Kasting 1988; 
Selsis et al. 2007b) because their contribution to a planet's albedo is expected to outweigh 
their contribution to the greenhouse effect. (A dense H2O atmosphere is already optically 
thick throughout most of the thermal-IR, so adding clouds has only a small effect on the 
outgoing IR radiation.) Conversely, CO2 ice clouds are expected to cause warming in a dense 
C02 atmosphere because they backscatter outgoing thermal-IR radiation more efficiently 
than they backscatter incoming visible/near-IR radiation (Forget & Pierrehumbert 1997). 
One can demonstrate the nature of these cloud influences using 1-D models, as was done 
in Selsis et al. (2007b). Making quantitative statements is difficult, however, because the 
warming or cooling effect of clouds depends on a host of parameters, including their heights, 
optical depths, particle sizes, and most importantly, fractional cloud coverage. Forget & 
Pierrehumbert (1997) obtained as much as 70 degrees of warming out of an optical depth 
10 CO2 cloud with 100% cloud cover, but that warming dropped by 30 degrees if fractional 
cloud cover was reduced to 75%. Realistic fractional cloud cover for condensation clouds is 
closer to 50%, because such clouds tend to form on updrafts, and approximately half the air 
in the troposphere is rising at any one time while the other half is descending. 

The best way to incorporate cloud effects in a climate calculation is to use a 3-D general 
circulation model (GCM). Attempts were made to explain warm early mars using such 3-D 
models (Forget et al. 2012) but none have yet succeeded. One can, however, do significantly 
better than in our 1-D model, and so further research in this area is warranted (Abe et al. 
2011; Wordsworth et al. 2011). 



Table 1: Habitable Zone distances around our Sun from our updated 1-D climate model. For 
comparison, estimates from Kasting et al. (1993) are also shown. 





Inner Habitable Zone 


Outer Habitable Zone 


Model 


Moist 
greenhouse 


Runaway 
greenhouse 


Recent Venus 


Maximum 
greenhouse 


Early Mars 


This paper 


0.99 AU 


0.97 AU 


0.75 AU 


1.70 AU 


1.77 AU 


Kasting et al. (1993) 


0.95 AU 


0.84 AU 


0.75 AU 


1.67 AU 


1.77 AU 



3.5. Habitable Zone Limits for Non-Earth-like planets 



In Table 2, we show the effect of surface gravities on the HZs of two planets. These 
planetary gravities were selected to encompass the mass range from Mars (gravity of 3.73 
ms~^) to a roughly 10 super-Earth (gravity of 25 ms~^). Both planets were assumed to 
have a 1 bar background N2 atmosphere. This may be unrealistic because proportionately 
more nitrogen is put on the smaller planet than the larger one; however, this allows direct 
comparison with Kasting et al. (1993). Table 2 shows that the habitability limits move 
slightly outward for a Mars-sized planet and inward for a super-Earth. This is because the 
column depth is larger for a Mars-sized planet, which increases the greenhouse effect (at the 
inner edge) and albedo (at the outer edge). Since the inner edge moves closer to the star 
for the super-Earth planet, while the outer edge changed little, we can conclude that, for a 
given surface pressure, larger planets have somewhat wider habitable zones than do small 
ones. 
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Table 2: Habitable Zones around our Sun for different planetary parameters. 



Inner Habitable Zone Outer Habitable Zone 


MnHpl 


Mnist 
greenhouse 
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1 OS^ ATI 


1 O'^S AU 


1 79 AU 


Earth 


0.99 AU 


0.97 AU 


1.70 AU 


Super-Earth** 


0.94 AU 


0.92 AU 


1.67 AU 


pCOs = 5.2 X 10-3 bart 


1.00 AU 


0.97 AU 




pCOs = 5.2 X 10-2 bar 


1.02 AU 


0.97 AU 




pCOs = 5.2 X 10-1 bar 


1.02 AU 


0.97 AU 




pCOa = 5.2 bar 


0.99 AU 


0.97 AU 





* Surface gravity = 3.73 m.s ^ 
** Surface gravity ~ 25 m.s"^ 

pC02 = 5.2 X 10^^ bar for our standard Earth model. Note that these CO2 pressures 

are not actual partial pressures; rather, they represent the surface pressure that would 

be produced if this amount of CO2 were placed in the atmosphere by itself. The 330 

ppmv of CO2 in our standard 1-bar atmosphere would produce a surface pressure of 

5.2 X 10^"* bar if the rest of the atmosphere was not present. When lighter gases such 

as N2 and O2 arc present, they increase the atmospheric scale height and cause CO2 

to diffuse upward, thereby lowering its partial pressure at the surface. 

We also performed sensitivity tests on the inner edge of the HZ by varying the amount 
of atmospheric CO2 (the outer edge calculation already factors in this change in CO2). It 
is quite possible that some terrestrial planets may have varying amount of CO2 because of 
different silicate weathering rates. As shown in Table. 2, changes in pC02 would not change 
the runaway greenhouse limit, as it is reached in an H20-dominated atmosphere. The moist 
greenhouse limit does change, as an increase in pC02 increases the surface temperature, and 
hence facilitates water loss. The maximum destabilization occurs at a PCO2 = 5.2 x 10-'^ 
bar approximately 10 times the present terrestrial PCO2 level (the critical distance, shown 
in bold in Table 2, is 1.00 AU). 
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This suggests that a 10-fold increase in CO2 concentration relative to today could push 
Earth into a moist greenhouse state (assuming a fully saturated atmosphere). By contrast, 
the maximum destabilization occurred at 1000 times the present CO2 level in Kasting et al. 
(1993). At larger pC02 values the increase in surface pressure outstrips the increase in the 
saturation vapor pressure of water, so the atmosphere becomes more stable against water 
loss (Kasting & Ackerman 1986). We conclude that planets with few tenths of a bar of PCO2 
have narrower HZs than planets like Earth on which pC02 is maintained at lower values by 
the carbonate-silicate cycle. 

4. Habitable Zones around Main-Sequence Stars 

The procedure described in the previous section to derive HZs around Sun can be used 
to estimate HZ boundaries around stars of different spectral types. A similar analysis was 
done by Kasting et al. (1993) for three stellar effective temperatures (7200 K, 5700 K and 
3700 K), which correspond to FO, GO, and MO spectral types. Selsis et al. (2007b) used a 
similar model to that of Kasting et al. (1993) and interpolated HZ distances to stars within 
this range of effective temperatures. Here, we compare our updated model results with these 
earlier studies and also extend the calculations to lower stellar effective temperatures to 
include M-dwarfs. Correctly calculating HZs of M-dwarfs is becoming increasingly impor- 
tant, as upcoming instruments such as Penn State's stabilized fiber-fed near-infrared (NIR) 
spectrograph Habitable Zone Planet Finder (HPF, Mahadevan et al. (2012)) and proposed 
missions such as Transiting Exoplanet Survey Satellite (TESS) will specifically search for 
low-mass planets around M-dwarfs. Furthermore, several rocky planets have already been 
found in the HZs of M-dwarfs (Bonfils et al. 2011; Vogt et al. 2012), and these objects may 
be good candidates for space-based characterization missions such as JWST. 

4.1. Habitable Zone Boundaries Around F, G, K and M Stars 

We considered stellar effective temperatures in the range 2600 K < Tg// < 7200 K, 
which encompasses F, G, K and M main-sequence spectral types. As input spectra for 
the HZ boundary calculations we used the "BT_Settl" grid of models^ (AUard et al. 2003, 
2007). These cover the needed wavelength range for climate models (0.23-4.54 /im), as well 
as the range of effective temperature (2600 K < T^fj < 70, 000 K) needed to simulate stellar 
spectra. Our comparison of the BT_Settl models with low-resolution IRTF data, and also 
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high-resolution CRIRES data on Barnard's star (from the CRIRES_POP hbrary^° (Lebzelter 
et al. 2012)), show that the models are quite good in reproducing the gross spectral features 
and energy distributions of stars, and will provide adequate input for our HZ calculations. 
For each star, the total energy flux over our climate model's spectral bands is normalized to 
1360 Wm~^ (the present solar constant for Earth) to simplify intercomparison. 

In Fig. 6 we compare the results of our inner and outer edge HZ model calculations 
for Sun to stars of different spectral types. Unless otherwise specified, we use the HITEMP 
2010 database for our inner edge calculations. The planetary albedo, shown in Fig. 6(a) 
(inner edge) and Fig. 6(c) (outer edge), of an Earth-like planet is higher if the host star is 
an F-star and lower if its primary is an M-star. The reason is that the Rayleigh scattering 
cross section (which is proportional to 1/A^) is on average higher for a planet around an 
F-star, as the star's Wien peak is bluer compared to the Sun. Second, H2O and CO2 have 
stronger absorption coefficients in the near-infrared than in the visible, so the amount of 
starlight absorbed by the planet's atmosphere increases as the radiation is redder (as is the 
case for an M-star). Both effects are more pronounced when the atmosphere is dense and 
full of gaseous absorbers. For a late M-star {T^ff = 2600 K) most of its radiation is peaked 
around 1 micron. Therefore, the minimal amount of Rayleigh scattering and the high near-IR 
absorption by the planet's atmosphere combine to generate extremely low planetary albedos. 

The changes in predicted planetary albedo can be translated into critical solar fluxes, 
as shown in Figs. 6(b) (inner edge) and 6(d) (outer edge). As discussed in §3.3, Seff goes 
through a minimum near the OHZ because the atmosphere becomes optically thick at all 
infrared wavelengths and, at the same time, the Rayleigh scattering due to CO2 condensation 
increases planetary albedo. Note that for a late M-star {T^ff = 2600 K) Rayleigh scattering 
never becomes an important factor, and hence Seff asymptotically reaches a constant value. 
The parameter Seff is directly calculated from our climate model and is dependent on 
the type of star considered. Therefore, we have derived relationships between HZ stellar 
fluxes (Seff) reaching the top of the atmosphere of an Earth-like planet and stellar effective 
temperatures (Teff) applicable in the range 2600 K < Teff < 7200 K. 

Seff = SeffQ + aT^ + hT^ + cTl + dTt (2) 

where = Teff — 5780 K and the coefficients are listed in Table 3 for various habitability 
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Fig. 6. — Habitable zone calculations from our climate model for stellar effective tempera- 
tures corresponding to F (T^// = 7200 K), G (Sun), K {Teff = 4800 K) and M (Te// = 3800 
K and 2600 K) spectral types. The inner edge results are shown in the top row (Figs. 6(a) 
& 6(b)) and the outer edge results are shown in the bottom row (Figs. 6(c) & 6(d)). 
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limits^^. The corresponding habitable zone distances can be calculated using the relation: 



where L/Lq is the luminosity of the star compared to the Sun. 

In Fig. 7 we compare HZ fluxes (and distances) calculated using Eqs.(2) & (3) for 
the moist greenhouse case, with Selsis et al. (2007b) 0% cloud results for different stellar 
effective temperatures. As shown in Fig. 7(a), for low Tg//, there are large differences 
at the inner edge (dashed and solid red curves) between the models. This is because the 
spectrum of low-mass stars shifts towards the longer wavelengths, resulting in more near-IR 
flux compared to high-mass stars. In both the models the atmosphere of a planet in the 
inner HZ is H20-dominated, and so there is strong absorption in the near-IR. Since our 
model uses the most recent HITEMP database which has more H2O lines in the near-IR, 
the moist greenhouse limit occurs at a lower flux (farther from the star). Also, Selsis et al. 
(2007b) assumed Tg// = 3700 K for stars with temperatures below this value. This amplifles 
the differences, as these low mass stars have their peak fluxes in near-IR. These differences 
in inner habitable zone boundaries may become important for present and upcoming planet 
flnding surveys around M-dwarfs such as MEARTH (Nutzman & Charbonneau 2008) and 
Penn State's HPF (Mahadevan et al. 2012), whose goal is to discover potentially habitable 
planets around M-dwarfs. 

The luminosity of a main sequence star evolves over time, and consequently the HZ 
distances (Eq.(3)) also change with time. One can calculate "continuous" HZ (GHZ) bound- 
aries within which a planet remains habitable for a specifled length of time (we chose 5 Gyr). 
In Fig. 7(b), we show GHZ boundaries as a function of stellar mass for both our model and 
Selsis et al. (2007b) model, taking into account the stellar evolutionary models of Baraffe et 
al. (1998) for solar metallicity stars. Noticeable differences between the two models are seen 
for low mass stars near the inner edge (as also seen in Fig. 7(a). The large differences in 
Seff from Fig. 7(a) do not appear as pronounced in Fig. 7(b) because it is a log scale and 
also because the GHZ distance is inversely proportional to the square root of S^ff (Eq.(3)). 

In order to assess the potential habitability of recently discovered exoplanets, equilibrium 
temperature {T^q) has been used as a metric (Borucki et al. 2011; Batalha et al. 2012). 
Assuming an emissivity of 0.9, the ranges of HZ boundaries are taken to be 185 K < Tgg < 303 
K (Kasting 2011b). We would like to stress that the stellar fluxes (5'e//) provide a better 



These coefficients can be downloaded in a machine readable format from the electronic version of the 
journal. A fortran code is also available to calculate HZ stellar fluxes. 
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metric for habitability than does Tgg. This is because Tgg involves an assumption about Ab 
(0.3, usually) that is generally not valid. This value of is good for present Earth around 
our Sun. For a planet around a late M-star, Ab can vary from 0.01 near the inner edge to 
0.1 at the outer edge (see Fig. 6), depending on its location. Similarly, Ab for an F-star 
can range in between 0.38 — 0.51 for the inner and outer edge, respectively. This changes 
the corresponding T^g, and so a uniform criterion for HZ boundaries based on Tf,g cannot be 
determined. 



Table 3: Coefficients to be used in Eq.(2) to calculate habitable stellar fluxes, and corresponding habitable 
zones (Eq.(3)), for stars with 2600 < T^ff < 7200 K. An ASCII file containing these coefficients can be 
downloaded in the electronic version of the paper. 
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2.9875 X 10-9 


1.5418 X 10-8 


1.7063 X 10-9 
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5. Discussion 

A straightforward application of the calculations presented in the previous sections is to 
apply them to currently known exoplanetary systems. Fig. 8 shows various habitable zone 
boundaries (expressed in terms of effective stellar flux) as a function of stellar Tg//. The or- 
bital parameters of the planets and stellar characteristics were obtained from exoplanets.org 
(Wright et al. 2011). The green-shaded habitable region is bounded by the moist greenhouse 
limit (inner edge) and the maximum greenhouse limit (outer edge). Several currently known 
terrestrial mass exoplanets that have been proposed to be in the HZ by various studies are 
also shown. An important insight that can be obtained from this flgure (and one that cannot 
be seen in the HZ distance plot. Fig. 7(b)) is that a terrestrial mass planet that lies within 
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Fig. 7. — Habitable zone fluxes (panel 7(a)) and corresponding distances (panel 7(b)) from 
our model (solid lines) compared to Selsis et al. (2007b) results (dashed lines) for different 
stellar effective temperatures. The inner HZ fluxes from our model (red solid lines) are 
for the moist greenhouse case and outer edge fluxes (blue solid lines) are for the maximum 
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the two vertical dashed-lines in the green shaded region is in the HZ irrespective of the type 
of star it is orbiting. The corresponding flux boundaries for which a rocky planet is 'defi- 
nitely' in the HZ are 0.842 and 0.42. Currently, two exoplanets fall within this region, HD 
40307g (Tuomi et al. 2012a) and Gl 581g (Vogt et al. 2010, 2012). The detection of planets 
and orbital parameters for G1581 is complicated by the low amplitudes of the signal, stellar 
activity, and possible red noise. We have included a '?' for Gl 581 system of planets in the 
plot, indicating that there is an active ongoing discussion in the literature about the number 
of planets in this system and their exact orbital parameters (Vogt et al. 2010; Pepe et al. 
2011b; Vogt et al. 2012; Baluev 2012). Furthermore, for stars with Teff < 5000 K, there 
is no clear distinction between runaway greenhouse and the moist-greenhouse limits. The 
reason is that for these stars, there are more photons available in the IR part of the spectrum, 
where H2O is a good absorber. Therefore a planet with H2O dominated atmosphere quickly 
goes into runaway once it reaches the moist-greenhouse limit. Note that another suggested 
HZ candidate planet, HD 85512b (Pepe et al. 2011a), receives more than 5 times the stellar 
flux received by our Earth, placing it even beyond the most liberal ('recent Venus') estimate 
of the inner edge. Hence, it is very likely that this planet is not in the HZ of its star. 

A question of importance to the exoplanet community is which HZ limits to choose 
when identifying potentially habitable planets. For current RV surveys and Kepler mission 
one should use the most conservative limits (moist greenhouse and maximum greenhouse), 
because this will give a lower limit on 77^, the fraction of sun-like stars that have at least 
one planet in the habitable zone (Luinine et al. 2008). If one is interested in designing a 
future flagship mission, such as Terrestrial Planet Finder (TPF) or Darwin, then using these 
conservative limits (which results in a lower limit on 77®) ensures that the telescope is not 
undersized. If, however, one was analyzing data obtained from such a telescope, the most 
optimistic limits (recent Venus and early Mars) should be used because one would not want 
to miss out on any potentially habitable planets. 

In Fig. 9, we show the incident stellar flux as a function of planetary mass for the 
currently known exoplanets. The masses are obtained from exoplanets.org when available. 
Also shown are habitable zone flux boundaries calculated from Eq.(2) for terrestrial mass 
planets (IM^ - lOM^). For the outer box (light grey), the upper bound on the flux is 
taken to be the moist greenhouse limit for a star with Tgjj = 7200 K and the lower bound 
is the maximum greenhouse limit for a star with Te// = 2600 K. These are indicated by 
the diagonally opposite points on the green shaded region in Fig. 8. For the inner box 
(dark grey), the flux limits are the dashed lines in the green shaded region of Fig. 8. The 
significance of this plot is that terrestrial planets in the inner box must be in the HZ, 
irrespective of the stellar spectral type. Mars, if it were more massive, would be in the HZ 
around any main sequence star with 2600 < T^ff < 7200 K. For planets that are outside 
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the dark grey region, but inside the hght grey (for example, Earth) one needs to know the 
host star's spectral type (or T^ff) to determine if that planet is in the HZ. Fig. 9 combines 
observable stellar and planetary parameters to further constrain HZ boundaries for extrasolar 
planets. 

Many of the currently known exoplanets have non-zero eccentricities, which can carry 
some of them (and their possible moons) in and out of the HZ. The incident stellar flux 
on these eccentric planets has extreme variations between periastron and apoastron ([(1 + 
e)/(l - e)Y). Williams & Pollard (2002) show that, provided that an ocean is present to 
act as a heat capacitor, it is primarily the time-averaged flux < SJjjr > that affects the 
habitability over an eccentric orbit (Kopparapu et al. 2009; Kopparapu & Barnes 2010). 
Mathematically: 



Here, Se/f is the effective flux from circular orbit (Eq.(2)). Planets with high orbital eccen- 
tricities (e ~> 0.1) have higher average orbital flux. This may help eccentric planets near 
the outer edge of the HZ maintain habitable conditions. However, obliquity variations can 
influence the geographical distribution of irradiation (Spiegel et al. 2008, 2009; Dressing et 
al. 2010) and may change habitable conditions. 

Earth itself appears to be perilously close to the moist greenhouse limit {Seff = 1.015, 
blue filled circle in Fig. 8). However, this apparent instability is deceptive, because the 
calculations do not take into account the likely increase in Earth's albedo that would be 
caused by water clouds on a warmer Earth. Furthermore, these calculations assume a fully 
saturated troposphere that maximizes the greenhouse effect. For both reasons, it is likely 
that the actual HZ inner edge is closer to the Sun than our moist greenhouse limit indicates. 
Note that the moist greenhouse in our model occurs at a surface temperature of 340 K. 
The current average surface temperature of the Earth is only 288 K. Even a modest (5-10 
degree) increase in the current surface temperature could have devastating affects on the 
habitability of Earth from a human standpoint. Consequently, though we identify the moist 
greenhouse limit as the inner edge of the habitable zone, habitable conditions for humans 
could disappear well before Earth reaches this limit. 

Additional uncertainty about habitability of planets around late-K and M stars. (Tg// <= 
4000 K) comes from the fact that planets within the HZs of these stars are expected to be 
tidally locked (Dole 1964; Peale 1977; Kasting et al. 1993; Dobrovolskis 2009). If the 
planet's orbital eccentricity is small, this can result in synchronous rotation, in which one 
side of a planet always faces the star (as the Moon does to the Earth). Climates of syn- 
chronously rotating planets are not well approximated by 1-D, globally averaged models. 
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Previous work has shown that such planets may indeed be habitable (Joshi et al. 1997; Joshi 
2003; Edson et al. 2003); however, systematic exploration of synchronously rotating planets 
in different parts of the HZ has not been attempted. Even before doing these calculations, 
we can predict that planets near the outer edge of the HZ, with their expected dense CO2 
atmospheres, should be more effective at transporting heat around to their night sides, and 
hence should have a better chance of being habitable. 

Given that survey like HPF and CARMENES will specifically target mid-late M dwarfs, 
our future work will include estimating the HZ boundaries of individual targets in detail. 
Ongoing work by our team (Terrien et al. 2012) is yielding low resolution NIR spectra 
from the infrared Telelescope Facility (IRTF) to be used to derive stellar metallicities as 
well as yielding more realistic flux distributions and temperatures for use in the modeling. 
We anticipate having this information for ~ 650 M dwarfs drawn from the J< 10 Lepine 
& Gaidos (2011)catalog, and have applied for time to observe ~ 300 more. Estimates of 
luminosities will be derived using photometric and spectroscopic distances for now (in cases 
where parallax measurements are absent), but eventually GAIA (Ferryman et al. 2001) will 
yield very precise parallaxes (and by extension precise luminosities) for all these target stars. 

Recent discoveries by both the Kepler mission and RV surveys have shown that planets 
can exist in stable orbits around multiple star systems (Doyle et al. 2011; Welsh et al. 2012; 
Orosz et al. 2012; Dumusque et al. 2012). The HZs of these stars could potentially host 
terrestrial planets, which are at the threshold of current detection techniques. Indeed, the 
discovery of Kepler 47c (Orosz et al. 2012) which is 4.6 times the size of the Earth's radii 
in the HZ is a step closer to discovering rocky planets in the HZ of multiple star systems. 
Dumusque et al. (2012) have recently published a possible detection of a 1.1 Earth mass 
(minimum) planet in a 3.236 day orbit around a Centauri system. So, this system should 
now be a prime target for further observations to discover habitable planets. Formation of 
dynamically stable terrestrial planets in the HZs of multiple star systems has been studied 
before (Whitmire et al. 1998; Holman & Wiegert 1999; Haghighipour & Raymond 2007) and 
several studies estimated HZ boundaries around these types of systems (Eggl et al. 2012a,b; 
Kane & Hinkel 2012) using Kasting et al. (1993) model. Our updated model results from 
Fig. 8 or Eqs.(2) and (3) could change these estimates significantly. 

Our new model results could also directly affect estimates of 77®. Recent analysis of 
Keplerdata. (Traub 2012) and RV surveys (Bonfils et al. 2011) concluded that ry® ~ 0.34—0.4. 
These values were based either on the Kasting et al. (1993) model (Traub 2012) or the 
Selsis et al. (2007b) results (Bonfils et al. 2011). Our new HZ limits could impact these 
estimates significantly. In particular, there are large differences between Selsis et al. (2007b) 
calculations and our model results for low mass stars. The estimate of 77© by Bonfils et al. 



(2011) is obtained by using Selsis et al. (2007b) relationships for planets orbiting M-stars. 
Thus, this value may need to be re-evaluated. 




Fig. 8. — Various cloud-free habitable zone (flux) boundaries for stars with different Tg//. 
The boundaries of the green-shaded region are determined by the moist-greenhouse (inner 
edge, higher flux values) & maximum greenhouse (outer edge, lower flux values). A planet 
that receives stellar flux bounded by the two dashed vertical lines is in the HZ irrespective 
of the stellar type. Some of the currently known exoplanets that are thought to be in the 
HZ by previous studies are also shown. The '?' for Gl 581 and Tau Ceti system of planets 
imply that there is an ongoing discussion about their existence. For stars with Teff ^ 5000 
K, there is no clear distinction between runaway and moist-greenhouse limit. 



6. Conclusions 

We have obtained new estimates for HZs around F, G, K and M main-sequence stars by 
(1) updating H2O and CO2 absorption coefficients in the Kasting et al. (1993) 1-D radiative- 
convective cloud-free climate model with the most recent LBL databases: HITRAN 2008 
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Fig. 9. — HZ boundaries combining observable stellar and planetary parameters. Terrestrial 
planets (O.3M0 - lOM^, y-axis limits) within the dark shaded region are in the HZ irrespective 
of the stellar type. For those planets in the light grey region, one needs to know the stellar 
type to determine if they are in the HZ. For example, the Earth would not be habitable if it 
received its current incoming amount of energy from a cooler-type star, but it is (obviously) 
habitable in its current orbit around a G-type star. For planets such as this, the star's energy 
distribution matters when considering habitability. 
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and HITEMP 2010, (2) correcting the H2O Rayleigh scattering coefficient, and (3) updating 
CO2 collision-induced absorption coefficients. These changes affect the inner and outer edges 
of the HZ, respectively. 

Our revised model predicts that the moist greenhouse limit for our Sun, which defines 
the inner edge of the HZ, is at 0.99 AU. The outer edge of the HZ, where gaseous CO2 
produces its maximum greenhouse effect, is at 1.70 AU. Although it appears that Earth is 
perilously close to the inner HZ edge, in reality, cloud feedback and low upper tropospheric 
relative humidity act to stabilize Earths climate. Theoretical studies and observational 
surveys that depend on these limits should use the updated values. We have also estimated 
HZ boundaries for M stars with T^ff as low as 2600 K, which are primary targets for ongoing 
surveys such as Habitable Zone Planet Finder and MEARTH to discover potential habitable 
planets. 

We also showed that the effective stellar flux provides a better criterion in determining 
the HZ limits, than equilibrium temperature. Accordingly we have derived a generalized 
expression to calculate these fluxes for stars of different spectral types. Our results show 
that some of the extrasolar planets that were previously thought to be within the HZ may 
not be in that region. Thus, our HZ estimates can be used to narrow the target list for even- 
tual characterization missions, such as JWST, to identify potential biomarkers on habitable 
planets. 
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